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99, 237 (1996), we introduced an abstract Newton-like equation on a general 
Lie algebra such that submanifolds fixed by the second-order Casimir operator 
are attracting set. The corresponding group parameters satisfy the nonlinear 
dynamical system having an attractor coinciding with the submanifold. In this 
work we discuss the case with the SU(2) x SU(2) group. The resulting second- 
order system in R 6 is demonstrated to exhibit chaotic behaviour. 
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One of the most important problems of the theory of dynamical systems relates 
to finding and clasifying attracting limit sets or simply attractors. Another 
important problem is how to study the system on its invariant manifold. In 
our previous work [1] we introduced a method for the construction of G- 
invariant nonlinear dynamical systems having an attracting set coinciding with 
a submanifold generated by the second-order Casimir operator referring to a 
given Lie group G. The example of the SU (2) group discussed therein led to 
the oscillatory dynamics on the orbit i.e. the sphere £2. As suggested in [1] 
such regular dynamics need not be the case for groups with higher-dimensional 
submanifolds. In this paper we examine the case with the £77(2) x SU(2) 
group. Namely, following the general scheme described in [1] we introduce the 
nonlinear second-order system in R 6 satisfied by the group parameters. We 
then show that the system exhibits chaotic behaviour on the sphere £ 5 which 
is a submanifold fixed by the second-order Casimir operator corresponding to 
the £77(2) x £77(2) group. 



2. THE NEWTON-LIKE EQUATION 



In this section we recall the abstract Newton-like equation on a general Lie 
algebra such that submanifolds fixed by the second-order Casimir operator 
are attracting set [1]. This equation generates the nonlinear dynamical system 
satisfied by the group parameters, having an attractor coinciding with the 
submanifold. Consider the following second-order differential equation on a 
Lie algebra g: 

pX + vX + pX + aY = e iX Ye~ iX , X(0) = X , X(0) = X , (2.1) 

where X(t): R — > Q is a curve in q, Y e Q is a fixed element, p, is, p, a e R 
and dot designates differentiation with respect to time. 

Let us assume that /1 ^ 0. Clearly, we can set without loose of generality 
p > 0. On demanding that (2.1) admits the solution on the submanifold fixed 
by the second-order Casimir operator of the form 

TrX 2 = const ^ 0, (2.2) 

rescaling t — > yfpt and setting = (5 we arrive at the following Newton-like 
equation: 

TrX 2 

X + /3X + -— X = e lX Ye~ tX -Y, X(0) = X , X(0) = X . (2.3) 
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TrX 2 = 2TrX X |(l - e~ pt ) + TrX 2 . (2.4) 

That is if (5 > then the solution to (2.1) approaches the sub manifold 

TrX 2 = f TrX X + TrX 2 . (2.5) 

It can be easily checked that for (3 < there is no solution on the submanifold. 
The only exception are the initial data such that 

TrX X = 0, (2.6) 

The more detailed analysis of eq. (2.3) and (2.4) is provided in [1]. Notice 
that in the general case of an n-dimensional compact Lie algebra the manifold 
defined by (2.2) is simply the sphere S n -i. 

3. DYNAMICS ON THE SU(2)xSU(2) GROUP 



Our aim now is to discuss the nonlinear dynamical system implied by the 
abstract Newton-like equation in the case of the SU (2) x SU (2) group. We 
first observe that due to the product structure of the considered group, the 
general elements of the Lie algebra X(t) and Y can be written in the form 

X(t) = X + (t)+X_(t), 

Y = Y + + Y_, (3.1a) 
where 

[X + , X_] = 0, [X + , Y_] = [X_,Y + ] = 0, 

TrX + X_ = TrX + X_ = 0. (3.1b) 
Notice that in view of (3.1) the manifolds (2.2) take the form 

Tr(X 2 +X 2 ) = const. (3.2) 
We now return to (2.3). From (2.3) and (3.1) we get 

*+ + ^ + W£m x + =e ' x+Y ^- Y+ ' 



X_ + PX^ + 772vX_ = e iX -Y_e~ iX - - YL, (3.3) 



Tr(X 2 +X_ 2 ' 
Tr(X 2 +X 2 ] 

X + (0) = X +0 , X_(0) = X_ , X + (0) = X +0 , X_(0) = X_ . 
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where cr = (<7i, a 2 , <7 3 ) and <jj, i — 1, 2, 3, are the Pauli matrices. Clearly, 
we can write the general elements of the Lie algebra X±(t) and Kj- of the Lie 
algebra as 



X + (t)=x + (t)-J + , X_(i)=x_(*).J_, 
F + = a + -J+, Yl = a -J_, 



(3.5) 



where x±(£): R — > R 3 , a-t is a constant vector of R 3 and the dot designates the 
inner product. Substituting (3.5) into (3.3) we obtain the following nonlinear 
system of second-order equations: 



x + + /3x + , 

X i 



X 



(cos |x + | — l)a + + 



x + /3x_ + 

xi 



xi + x! 



x + (0) 



H-0? 



2 x = (cos |x_| — l)a_ + 

x_o, x+(0) = x +0 



X 

x_(0) 



sin 


x +l 


X 

sin 


+ 

X 






X 







a + x x + + (1 — cos |x + |) 



v a + -x + jx^ 

A 
(a_-x_)x_ 



a x x + (1 — cos |x_|)- 

xl 

x_(0) = x_ , 



where a ± x x± designates the vector product of vectors a ± and x±, and |x±| = 
y^x4 is the norm of the vector x±. 

Notice that the manifolds given by (3.2) and (3.5) are the five- dimensional 
spheres 



x + + x_ = const. 



(3.7) 



(3.6) 



We now examine the asymptotic behaviour of the system (3.6). First we ob- 
serve that (3.6) implies 



<9x 



x+ + tt— -x_ = -Q{3 - 2-^ 



X -X 



dk. 



(3.8) 



Thus as the solution to (3.6) when (3 > approaches the sphere (3.7), the 
system (3.6) becomes dissipative one with exponential contraction of a vol- 
ume element. Furthermore, it can be easily checked that the £77(2) x SU(2) 
realization of (2.4) takes the form 

x^ + x 2 = 2(x +0 -x +0 + x_ -x_ )i(l - e~ pt ) + x% + x 2 . (3.9) 

Thus whenever the initial data satisfy the inequality 



|(x +0 -x +0 + x_ -x_o) + x 2 +0 + x 2 > 0, 



(3.10) 
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x; + xl = |(x + o-x+o + x_o-x_o) + Ko + x -o- I 3 - 11 ) 

Evidently, the sphere S5 is an invariant set for the initial conditions such that 

x +0 -x +0 + x„ -x_o = (3.12) 

and arbitrary (3. The solutions corresponding to the remaining initial data do 
not approach the sphere S 5 . Indeed, if (3 > and 

|(x +0 -x +0 + x_o-x_o) +x+ + x 2 = 0, x + o-x + o + x_o-x_ 7^ 0, (3.13) 

then the solutions to (3.6) tend asymptotically to the singular point x + = 
0, x_ = 0. Furthermore, if (3 > and 

|(x +0 -x +0 + x_ -x_o) + x% + x 2 < 0, (3.14) 

then the singular point is approached after a finite period of time 

(3 \ 2(x +0 -x +0 + x_o-x_+o)/ 
Finally, for (3 < and x +0 -x +0 + x_ -x_o 7^ 0, the trajectories go to infinity. 

We now return to (3.6). Notice that in view of the covariance of (3.3) with 
respect to the group transformations we can set in (3.6) without loose of 
generality 

a + = (0, 0, acosa), a + = (0, 0, asina). (3.16) 
Thus we finally arrive at the system 

x^ + x 2 sin|x±| ,,a± 3 x ±1 x± 3 

X±l + PX±! + -5— T X ±1 = : 7—a± 3 X± 2 + (1 - COS X± ) 2 , 

X+ + XI |X ± | X+ 

x^ + x 2 sin|x±| l ,a ±3 x ±2 x ±3 

x± 2 + px± 2 + — T x±2 = — ^a ±3 x ± i + (1 - cos x± ) - 2 , 

x+ + xi |x±| x| 

x ±3 + /*c ±3 + — 2-a; ± 3 = (cos |x±| - l)a ±3 + (1 - cos |x ± |) - 2 — , 

X_|_ + X X_|_ 

x±(0)=x ±0 , x±(0)=x ±0 , (3.17) 

where a +3 = a cos a and a_3 = a sin a. Notice that the parameter a can be 
restricted to the interval [0, 7r/2]. Indeed, taking into account (3.16) we find 
that the transformation a — > a + 7r/2 leads to 

a +3 — > -a_ 3 , a_ 3 — > a +3 , (3.18) 

or in view of (3.17) 

x — > x+, x+i — > — x+2 — > x_2, x +3 — > -x_ 3 . (3.19) 
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should also be noted that by virtue of the following transformation law of the 
equation (2.3) under the scaling t — > At: 

TrX 2 

X + \{3X + T^r^X = e iX X 2 Ye- iX - X 2 Y, (3.20) 

we have actually two bifurcation parameters: /3 2 /a, where a 7^ 0, and a. 

In order to discuss the symmetries of (3.17) consider the original abstract 
equation (2.3). Evidently, (2.3) is invariant under the transformations referring 
to the stability group of Y, that is transformations leaving Y unchanged. 
Hence, making use of (3.3) and (3.5) we find that solutions to (3.17) have the 
form invariant under rotations about x +3 and rr_ 3 axes such that 

x' ±1 = cosip±x±i + sirnp±x±2, 
x' ±2 = — sin (f±x±i + cosip±x± 2 , 

x'± 3 = x ±3 , ^ ± G[0,2vr). (3.21) 

This suggests that whenever the initial condition x +0 , x_ , x +0 and x_ corre- 
sponds to an attractor A then the basin of attraction of A contains the points 
x' +0 , x'_ , x' +0 and x'_ of the form 

x'± 10 = cos(p±x± w + sinv? ± a; ± 20, 
x '±20 = - sin ¥± x ±w + cosy? ± x ±2 o, 

^±30 = x ±30, (3.22a) 
x' ±w = cos (p±x± 10 + siny2 ± x ± 20, 
^±20 = ~ sin(p ± x± w + cosv5±x ±2 o, 

i±3o = ^±30, <P± e [0,2tt). (3.22b) 

Thus the basin of attraction contains two circles given by (3.22a) such that 

x +0 = x +0 , x _ = x_ , (3.23) 

and two vector fields (3.22b) obtained from x±o by rotating this about x +3 
and X-3 axis, respectively, by the angle referring to the position of the corre- 
sponding point of the circle. 

In figures 1 and 2 we show examples of strange attractors from numerical 
integration of the system (3.17). As expected these attractors are symmetric 
under rotations about x +3 and x_ 3 axes. 

Fig. 1, 2 
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Fig. 1 there exists a quasiperiodic trajectory. 



Fig. 3 



A look at Fig. 3 is enough to conclude that in the case with the attractor 
from Fig. 1 we deal with the quasiperiodicity to chaos transition like in the 
Ruelle-Takens-Newhouse scenario [2,3]. In the case of the attractor from Fig. 
2 we have most probably the new scenario. 



Fig. 4-6 



Namely, it turns out that there exist in the parameter space of the system 
(3.17) two nearby quasiperiodic orbits shown in Fig. 4. The computer simula- 
tions suggest that these quasiperiodic orbits are separated by an infinitesimal 
perturbation of the bifurcation parameter (3. As a consequence of the infinites- 
imal nature of the perturbation the irregular transitions occur between the 
chaotic attractors arising from the perturbation of the quasiperiodic orbits 
from Fig. 4a and Fig. 4b, respectively. These attractors are shown separately 
in Fig. 6. We remark that the attractors have the same form as those shown 
in Fig. 5 arising in the transient chaos before reaching the quasiperiodic state. 
The authors did not find such scenario of transitions between chaotic attrac- 
tors in the literature. As the bifurcation parameter decreases and approaches 
the value corresponding to the attractor from Fig. 2 the frequency of transi- 
tions between the two attractors increases. On the other hand, the decay of the 
bifurcation parameter leads to the deformation of the attractors. As a result 
of these two combined processes we arrive at the attractor illustrated in Fig. 
2. In this sense that attractor arise from the perturbation of two quasiperiodic 
states. 

4. CONCLUSION 



In the present work we have studied the particular SU (2) x SU{2) realization 
of the abstract Newton-like equation having attractors coinciding with the 
submanifolds fixed by the second-order Casimir operator referring to a given 
Lie group. The resulting second-order system in R 6 has been shown to exhibit 
chaotic behaviour. The concrete examples of the chaotic attractors with the 
SU(2) x SU(2) symmetry have been provided. One of them (see Fig. 2 and 
discussion below) is related to the most probably new scenario of transitions 
between chaotic attractors and effectively forming one by perturbation of two 
quasiperiodic states. We point out that the case with the SU (2) x SU (2) real- 
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gated system (3.17) is the first example of the chaotic dynamical system on the 
manifold determined by the continuous Lie-group structure. We remark that 
the importance of such an example was indicated in [4-6]. Furthermore, the 
knowledge of chaotic attractors with fixed symmetry and known the invariant 
measure like those introduced in this work would be useful for testing methods 
of detecting symmetry of chaotic attractors such as for example the method of 
detectives [7] . Whenever the theory of groups appears to lead to some insight 
into the nature of chaos then it seems that the observations introduced herein 
would be an important point of departure in solving numerous problems. 
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Fig. 1. The system (3.17) with /3 = 0.1, a = 0.5, a = 0.2, x +0 = (1, 1, 1.5), 
x_ = (1,1,1.5), x + o = (-0.1,0.1,0.1) and x_ = (0.1,0.1,0.1). Left: the 
projection of the Poincare section of the attractor on the (x + i,x + s) plane. 
This section and the following ones in the case with the projection on the 
(x + i,x +3 ) plane are defined by the hyperplane n(x — x^) = 0, where n is 
the normal vector and x^ is a point of the hyperplane, n, x^°\ x £ R 12 , 
n = (0, 1, 0, . . . , 0) and x^ = 0. Right: the projection of the Poincare section 
of the attractor on the (x-i,X-s) plane. This section and the following ones 
in the case with the projection on the (x-i,xs) plane are defined by the 
hyperplane n(x - = 0, where n = (0, 0, 0, 0, 1, 0, . . . , 0) and x^ = 0. 
The Lyapunov exponents are: Ai = 0.02, A 2 = A 3 = A 4 = A 5 = ±0.00, 
A 6 = -0.01, A 7 = -0.06, A 8 = -0.08, A 9 = -0.09, Ai = -0.10, An = -0.12 
and A12 = —0.16. The Lyapunov dimension is 6.1. 

Fig. 2. The system (3.17) with (5 = 0.1, a = 1, a = tt/4, x +0 = (2, 2, 2), x_ = 
(2,2,2), x +0 = (0.1,0.1,-1) and x_ = (1,1,1). Left: the projection of the 
Poincare section of the attractor on the (x+i, x + %) plane. Right: the projection 
of the Poincare section of the attractor on the £-3) plane. The Lyapunov 
exponents are: Ai = 0.06, A 2 = 0.01, A 3 = A 4 = A 5 = A 6 = ±0.00, A 7 = -0.04, 
A 8 = -0.07, A 9 = -0.10, A 10 = -0.11, An = -0.14 and A 12 = -0.21. The 
Lyapunov dimension is 7.4. 

Fig. 3. The system (3.17) with f3 = 0.117. The remaining parameters and 
the initial condition are the same as in Fig. 1. Left: the projection of the 
Poincare section of the quasiperiodic attractor on the (x + i,x + z) plane. Right: 
the projection of the Poincare section of the quasiperiodic attractor on the 
(x-i, x-s) plane. 

Fig. 4. The system (3.17) with a) (5 = 0.1102766 and b) (5 = 0.1102767. 
The remaining parameters and the initial condition are the same as in Fig. 2. 
The projection of the Poincare section of the quasiperiodic attractors on the 
(x +1 ,x +3 ) plane. 

Fig. 5. The transient chaos in the system (3.17) with a) (3 — 0.1102766 and 
b) P — 0.1102767 before reaching the quasiperiodic states from Fig. 4. The 
remaining parameters and the initial condition are the same as in Fig. 2. The 
projection of the Poincare section of the attractors on the (x + i,x + s) plane. 

Fig. 6. The system (3.17) with f3 = 0.1099. The remaining parameters and 
the initial condition are the same as in Fig. 2. For such data the irregular 
transitions appear between the attractor from the left and the attractor from 
the right. In both figures the projection is shown of the Poincare section of 
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